A DISCRETE THEORY OF CONNECTIONS ON PRINCIPAL BUNDLES 
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Abstract. Connections on principal bundles play a fundamental role in expressing the equations of motion 
for mechanical systems with symmetry in an intrinsic fashion. A discrete theory of connections on principal 
bundles is constructed by introducing the discrete analogue of the Atiyah sequence, with a connection 
corresponding to the choice of a splitting of the short exact sequence. Equivalent representations of a discrete 
connection are considered, and an extension of the pair groupoid composition, that takes into account the 
principal bundle structure, is introduced. Computational issues, such as the order of approximation, are 
also addressed. Discrete connections provide an intrinsic method for introducing coordinates on the reduced 
space for discrete mechanics, and provide the necessary discrete geometry to introduce more general discrete 
symmetry reduction. In addition, discrete analogues of the Levi-Civita connection, and its curvature, are 
introduced by using the machinery of discrete exterior calculus, and discrete connections. 
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1. Introduction 

One of the major goals of geometric mechanics is the study of symmetry, and its consequences. An 
important tool in this regard is the non-singular reduction of mechanical systems under the action of free 
and proper symmetries, which is naturally formulated in the setting of principal bundles. 

The reduction procedure results in the decomposition of the equations of motion into terms involving the 
shape and group variables, and the coupling between these are represented in terms of a connection on the 
principal bundle. 

Connections and their associated curvature play an important role in the phenomena of geometric phases. 
A discussion of the history of geometric phases can be found in Berry [1990]. Shapere and Wilczek [1989] 
is a collection of papers on the theory and application of geometric phases to physics. In the rest of this 
section, we will survey some of the applications of geometric phases and connections to geometric mechanics 
and control, some of which were drawn from Marsden [1994, 1997], Marsden and Ratiu [1999]. 
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The simulation of these phenomena requires the construction of a discrete notion of connections on 
principal bundles that is compatible with the approach of discrete variational mechanics, and it towards this 
end that this chapter is dedicated. 

Falling Cat. Geometric phases arise in nature, and perhaps the most striking example of this is the falling 
cat, which is able to reorient itself by 180°, while remaining at zero angular momentum, as show in Figure 1. 

The key to reconciling this with the constancy of the angular momentum is that angular momentum 
depends on the moment of inertia, which in turn depends on the shape of the cat. When the cat changes it 
shape by curling up and twisting, its moment of inertia changes, which is in turn compensated by its overall 
orientation changing to maintain the zero angular momentum condition. The zero angular momentum 
condition induces a connection on the principal bundle, and the curvature of this connection is what allows 
the cat to reorient itself. 

A similar experiment can be tried on Earth, as described on page 10 of Vedral [2003]. This involves 
standing on a swivel chair, lifting your arms, and rotating them over your head, which will result in the chair 
swivelling around slowly. 

Holonomy. The sense in which curvature is related to geometric phases is most clearly illustrated by 
considering the parallel transport of a vector around a curve on the sphere, as shown in Figure 2. 

Think of the point on the sphere as representing the shape of the cat, and the vector as representing its 
orientation. The fact that the vector experiences a phase shift when parallel transported around the sphere 
is an example of holonomy. In general, holonomy refers to a situation in geometry wherein an orthonormal 
frame that is parallel transported around a closed loop, back to its original position, is rotated with respect 
to its original orientation. 

Curvature of a space is critically related to the presence of holonomy. Indeed, curvature should be thought 
of as being an infinitesimal version of holonomy, and this interpretation will resurface when considering the 
discrete analogue of curvature in the context of a discrete exterior calculus. 

Foucault Pendulum. Another example relating geometric phases and holonomy is that of the Foucault 
pendulum. As the Earth rotates about the Sun, the Foucault pendulum exhibits a phase shift of A9 = 
27T cos (x (where a is the co- latitude) . This phase shift is geometric in nature, and is a consequence of 
holonomy. If one parallel transports an orthonormal frame around the line of constant latitude, it exhibits 
a phase shift that is identical to that of the Foucault pendulum, as illustrated in Figure 3. 

True Polar Wander. A particular striking example of the consequences of geometric phases and the 
conservation of angular momentum is the phenomena of true polar wander, that was studied by Goldreich 
and Toomre [1969], and more recently by Leok [1998]. It is thought that some 500 to 600 million years ago, 
during the Vendian-Cambrian transition, the Earth, over a 15-million-year period, experienced an inertial 
interchange true polar wander event. This occurred when the intermediate and maximum moments of inertia 
crossed due to the redistribution of mass anomalies, associated with continental drift and mantle convection, 
thereby causing a catastrophic shift in the axis of rotation. 

This phenomena is illustrated in Figure 4, wherein the places corresponding to the North and South poles 
of the Earth migrate towards the equator as the axis of rotation changes. 

Geometric Control Theory. Geometric phases also have interesting applications and consequences in 
geometric control theory, and allow, for example, astronauts in free space to reorient themselves by changing 
their shape. By holding one of their legs straight, swivelling at the hip, and moving their foot in a circle, 
they are able to change their orientation. Since the reorientation only occurs as the shape is being changed, 
this allows the reorientation to be done with extremely high precision. Such ideas have been applied to the 
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Figure 3: Geometric phase of the Foucault pendulum. 



control of robots and spacecrafts; see, for example, Walsh and Sastry [1993]. The role of connections in 
geometric control is also addressed in-depth in Marsdcn [1994, 1997]. 

One of the theoretical underpinnings of the application of geometric phases to geometric control was 
developed in Montgomery [1991] and Marsden et al. [1990], in the form of the rigid-body phase formula, 

A0 = Ai/ "i> + 2H v T ) = - A + nnr> 
w\\ Ud J \w\ 

the geometry of which is illustrated in Figure 5. 

An example that has been studied extensively is that of the satellite with internal rotors, with a configu- 
ration space given by Q = SE(3) x S 1 x S 1 x S 1 , and illustrated in Figure 6. 

The generalization of the rigid-body phase formula in the presence of feedback control is particularly 
useful in the study and design of attitude control algorithms. 
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Figure 4: True Polar Wander. Red axis corresponds to the original rotational axis, and the gold axis 
corresponds to the instantaneous rotational axis. 




Figure 5: Geometry of rigid-body phase. 
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rigid carrier 




spinning rotors 
Figure 6: Rigid body with internal rotors. 

2. General Theory of Bundles 

Before considering the discrete analogue of connections on principal bundles, we will review some basic 
material on the general theory of bundles, fiber bundles, and principal fiber bundles. A more in-depth 
discussion of fiber bundles can be found in Steenrod [1951] and Kobayashi and Nomizu [1963]. 

A bundle Q consists of a triple (Q, S, tt), where Q and S are topological spaces, respectively referred to as 
the bundle space and the base space, and tt : Q — ► S is a continuous map called the projection. We may 
assume, without loss of generality, that tt is surjective, by considering the bundle over the image tt(Q) C 5. 

The fiber over the point x 6 S, denoted F x , is given by, F x = tt^ 1 (x). In most situations of practical 
interest, the fiber at every point is homeomorphic to a common space F, in which case, F is the fiber of the 
bundle, and the bundle is a fiber bundle. The geometry of a fiber bundle is illustrated in Figure 7. 



Q 




~~ s 



Figure 7: Geometry of a fiber bundle. 

A bundle (Q, S, tt) is a G -bundle if G acts on Q by left translation, and it is isomorphic to {Q,Q/G,ttq/g), 
where Q/G is the orbit space of the G action on Q, and ttq/g is the natural projection. 

If G acts freely on Q, then (Q, S, tt) is called a principal G-bundle, or principal bundle, and G is its 
structure group. G acting freely on Q implies that each orbit is homeomorphic to G, and therefore, Q is 
a fiber bundle with fiber G. 

To make the setting for the rest of this chapter more precise, we will adopt the following definition of a 
principal bundle, 
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Definition 2.1. A principal bundle is a manifold Q with a free left action, p : G x Q —> Q, of a Lie group 
G, such that the natural projection, 77 : Q — > Q/G, is a submersion. The base space Q/G is often referred to 
as the shape space S, which is a terminology originating from reduction theory. 

We will now consider a few standard techniques for combining bundles together to form new bundles. 
These methods include the fiber product, Whitney sum, and the associated bundle construction. 

Fiber Product. Given two bundles with the same base space, we can construct a new bundle, referred to 
as the fiber product, which has the same base space, and a fiber which is the direct product of the fibers 
of the original two bundles. More formally, we have, 

Definition 2.2. Given two bundles 77 j : Qi — ► S, i = 1, 2, the fiber product is the bundle, 

Tl X S 7T 2 : Qi x s Qi -> S, 

where Q\ Xg Q 2 is the set of all elements (91,92) E Qi x Q 2 such that 771(91) = 772(92); and the projection 
77i X5772 is naturally defined by tt\ Xg 772(91, 92) = 771(91) = 7r 2 (g 2 ). The fiber is given by (771 Xq 7r 2 ) _1 (a;) = 
Trf 1 ^) x ir^ix). 

Whitney Sum. The Whitney sum combines two vector bundles using the fiber product construction. 

Definition 2.3. Given two vector bundles n : Vi — ► Q, i = 1, 2, with the same base, their Whitney sum 
is their fiber product, and it is a vector bundle over Q, and is denoted V\ © V 2 . This bundle is obtained by 
taking the fiberwise direct sum of the fibers of V\ and V 2 • 

Associated Bundle. Given a principal bundle, 77 : Q — > Q/G, and a left action, p : G x M — > M, of the 
Lie group G on a manifold M, we can construct the associated bundle. 

Definition 2.4. An associated bundle M with standard fiber M is, 

M = Qx G M = (QxM)/G, 

where the action of G on Q x M is given by g(q,m) — (gq,gm). The class (or orbit) of {q,m) is denoted 
[q,m]G or simply [q, m]. The projection ttm '■ Q x G M — > Q/G is given by, 

■km : {[q,m] G ) = n(q) , 
and it is easy to check that it is well-defined and is a surjective submersion. 

3. Connections and Bundles 

Before formally introducing the precise definition of a connection, we will attempt to develop some intuition 
and motivation for the concept. As alluded to in the introduction to this chapter, a connection describes the 
curvature of a space. In the classical Riemannian setting used by Einstein in his theory of general relativity, 
the curvature of the space is constructed out of the connection, in terms of the Christoffcl symbols that 
encode the connection in coordinates. 

In the context of principal bundles, the connection provides a means of decomposing the tangent space 
to the bundle into complementary spaces, as show in Figure 8. Directions in the bundle that project to 
zero on the base space are called vertical directions, and a connection specifies a set of directions, called 
horizontal directions, at each point, which complements the space of vertical directions. 
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Figure 8: Geometric phase and connections. 



In the rest of this section, we will formally define connections on principal bundles, and in the next section, 
discrete connections will be introduced in a parallel fashion. 

Short Exact Sequence. This decomposition of the tangent space TQ into horizontal and vertical subspaces 
yields the following short exact sequence of vector bundles over Q, 



0- 



^VQ 



-> TQ ■ 



-+0, 



where VQ is the vertical subspace of TQ, and ir*TS is the pull-back of TS by the projection it : Q — > S. 

Atiyah Sequence. When the short exact sequence above is quotiented modulo G, we obtain an exact 
sequence of vector bundles over S, 







-4 TQ/G — ^ TS ► , 



which is called the Atiyah sequence (see, for example Atiyah [1957], Almeida and Molino [1985], Mackenzie 
[1995]). Here, g is the adjoint bundle, which is a special case of an associated bundle (see Definition 2.4). 
In particular, 

a = Q xg0 = (Q x q)/g, 

where the action of G on Q x g is given by g(q, £) = (gq, Ad g £), and 7r : g — ► S is given by 7r fl ([g, £\g) = ^(q)- 
The maps in the Atiyah sequence, i : (Q x Q)/G — > TQ/G and 7r» : TQ/G — ► TS, are given by 

i{[Q,$G) = Ko(«)]c, 

and 

tt*(K]g) = r7r(u 9 ). 

Connection 1-form. Given a connection on a principal fiber bundle tt : Q — > Q/G, we can represent this as 
a Lie algebra- valued connection 1-form, A : TQ — ► g, constructed as follows (see, for example, Kobayashi 
and Nomizu [1963]). Given an element of the Lie algebra (eg, the infinitesimal generator map £ h- > £q 
yields a linear isomorphism between g and V^Q for each q E Q. For each w g G T«jQ, we define ^l(w g ) to be 
the unique £ G g such that £q is equal to the vertical component of v q . 

Proposition 3.1. The connection 1-form, A : TQ — > g, o/ a connection satisfies the following conditions. 
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(1) The 1-form is G-equivariant, that is, 



A o TL g = Ad,, oA , 



for every g G G, where Ad denotes the adjoint representation of G in g. 
(2) The 1-form induces a splitting of the Atiyah sequence, that is, 

for every £ G g. 

Conversely, given a g-valued 1-form A on Q satisfying conditions 1 and 2, there is a unique connection in 
Q whose connection 1-form is A. 

Proof. See page 64 of Kobayashi and Nomizu [1963]. □ 

Horizontal Lift. The horizontal lift of a vector field X G X(S) is the unique vector field X h G X(Q) 
which is horizontal and which projects onto X, that is, Tir q (Xq) = X v ^ for all q G Q. The horizontal lift 
is in one-to-one correspondence with the choice of a connection on Q, as the following proposition states. 

Proposition 3.2. Given a connection in Q, and a vector field X G %(S), there is a unique horizontal lift 
X h of X . The lift X h is left-invariant under the action of G. Conversely, every horizontal vector field X h 
on Q that is left-invariant by G is the lift of a vector field X G X(S). 

Proof. See page 65 of Kobayashi and Nomizu [1963]. □ 
Connection as a Splitting of the Atiyah Sequence. Consider the continuous Atiyah sequence, 

>gr^TQ/G^^TS >0 

(tti,A) X h 

We see that the connection 1-form, A : TQ — > g, induces a splitting of the continuous Atiyah sequence, since 

(iri,A)oi([q,^] G ) = (tti, AM Q (q)] g ) = [q,A{t Q (q))} G = [q,Z] G , for all q G Q,£ G g, 

which is to say that (tti, A) o i = lg. Conversely, given a splitting of the continuous Atiyah sequence, we can 
extend the map, by equivariance, to yield a connection 1-form. 

The horizontal lift also induces a splitting on the continuous Atiyah sequence, since, by definition, the 
horizontal lift of a vector field X G X(S) projects onto X, which is to say that 7r* oX h = Its- The horizontal 
lift and the connection are related by the fact that 

!tq/g = * ° {m,A) + X h o 71-*, 

which is a simple consequence of the fact that the two splittings arc part of the following commutative 
diagram, 

► g TQ/G TS ► 



X" 

7T2 



Its 



> g 8®TS TS > 



7Tl 



12 



where aj. is an isomorphism. The isomorphism is given in the following lemma. 
Lemma 3.3. The map a A : TQ /G —> g © TS defined by 

(*A([q, Q\g) = [q, A(q, q)] G © Tir{q, q), 
is a well-defined vector bundle isomorphism. The inverse of a A is given by 

oi 1 (Mcffi(it.i)) - [(x,x)$ + Zq]g- 
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Proof. See page 15 of Cendra et al. [2001]. □ 

This lemma, and its higher-order generalization, that identifies T^Q/G with T^S Xs 2g, is critical in 
allowing us to construct the Lagrange-Poincare operator, which is an intrinsic method of expressing the 
reduced equations arising from Lagrangian reduction. 

In the next section, we will develop the theory of discrete connections on principal bundles in a parallel 
fashion to the way we introduced continuous connections. 

4. Discrete Connections 

Discrete variational mechanics is based on the idea of approximating the tangent bundle TQ of Lagrangian 
mechanics with the pair groupoid Q x Q. As such, the purpose of a discrete connection is to decompose the 
subset of Q x Q that projects to a neighborhood of the diagonal oi S x S into horizontal and vertical spaces. 

The reason why we emphasize that the construction is only valid for the subset of Q x Q that projects 
to a neighborhood of the diagonal of S x S is that there are topological obstructions to globalizing the 
construction to all of Q x Q except in the case that Q is a trivial bundle. 

One of the challenges of dealing with the discrete space modelled by the pair groupoid Q x Q is that it 
is not a linear space, in contrast to TQ. As we shall see, the standard pair groupoid composition is not 
sufficient to make sense of the notion of an element (qo, q\) G Q x Q being the composition of a horizontal 
and a vertical clement. We will propose a natural notion of composing an clement with a vertical element 
that makes sense of the horizontal and vertical decomposition. 

In the subsequent sections, we will use the discrete connection to extend the pair groupoid composition 
even further, and explore its applications to the notion of curvature in discrete geometry. 

Intrinsic Representation of the Tangent Bundle. The intuition underlying our construction of discrete 
horizontal and vertical spaces is best developed by considering the intrinsic representation of the tangent 
bundle. This representation is obtained by identifying a tangent vector at a point on the manifold with the 
equivalence class of curves on the manifold going through the point, such that the tangent to the curve at 
the point is given by the tangent vector. This notion is illustrated in Figure 9. 




Figure 9: Intrinsic representation of the tangent bundle. 



Given a vector v q G TQ, we identify it with the family of curves q : M — > Q, such that q(0) — q, and 
q(0) = v. The equivalence class [•] identifies curves with the same basepoint, and the same velocity at the 
basepoint. 

With this representation, it is natural to consider (q 0} qi) G Q x Q to be an approximation of [(?(•)] 
TQ, in the sense that, 

<?o = q(0), qi = q(h), 
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for some fixed time step h, and where q(-) is a representative curve corresponding to v q in the intrinsic 
representation of the tangent bundle. 

4.1. Horizontal and Vertical Subspaces of Q x Q. Recall that the vertical subspace at a point q, denoted 
V q , is given by 

V q = {v q e TQ | = 0} = {£ Q I £ G 0}. 

Notice that the vertical space is precisely that subspace of TQ which maps under the lifted projection map 
to the embedded copy of S in TS. We proceed in an analogous fashion to define a discrete vertical subspace 
at a point q. 

The natural discrete analogue of the lifted projection map 7r* is the diagonal action of the projection map 
on Q x Q, (ir, 7r) : Q x Q — > Q x Q, where (qo, qi) >— > (7rg , nqi)- This is because 

7r.K) = 7r.([«(-)]) = [7r(g(-))]. 
In the same way that we embed 5 into TS* by the map x t-^> [x] = X , S naturally embeds itself into the 
diagonal of S x S, x i— > (a;, a;) = esxs, which we recall is the identity subspace of the pair groupoid. 

The alternative description of the vertical space is in terms of the embedding of Q x g into TQ, by 
(q, £) i — ► £q(<z), using the infinitesimal generator construction, 

= [exp(^)g]. 

In an analogous fashion, we construct a discrete generator map, which is given in the following definition. 
Definition 4.1. The discrete generator is the map i: QxG^QxQ, given by 

i(q,g) = (q,gq), 

which we also denote by i q {g) = i(q,g) = (q,gq). 

Then, we have the following definition of the discrete vertical space. 
Definition 4.2. The discrete vertical space is given by 

Ver q = {{q,q') eQxQ\ {n,n)(q,tf) = e SxS } 
= {i q {g) \geG}. 

This is the discrete analogue of the statement Ver g = {v q £ TQ \ ■n*{v q ) = 0} = {£q | ^ E g}. 

Since the pair groupoid composition is only defined on the space of composable pairs, we need to extend 
the composition to make sense of how the discrete horizontal space is complementary to the discrete vertical 
space. In particular, we define the composition of a vertical element with an arbitrary element of Q x Q as 
follows. 

Definition 4.3. The composition of an arbitrary element (qo,qi) £ Q x Q with a vertical element is given 
by 

iqo(g) ■ (<7o,<7i) = (e,g)(q ,qi) = (qo,gqi). 

An elementary consequence of this definition is that it makes the discrete generator map a homomorphism. 

Lemma 4.1. The discrete generator, i q , is a homomorphism. This is a discrete analogue of the statement 
in the continuous theory that (£ + x)q = £q + XQ- 

Proof. We compute, 

iq(g) ■ iq( h ) = l q{g) ' (?) h q) 
= (e,g)(q,hq) 

= (q,ghq) 

= iq{gh)- 
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Therefore, i q is a homomorphism. □ 

If we define the G action on Q x G to be h(q, g) = (hq, hgh -1 ), we find that the composition of a vertical 
element with an arbitrary element is G-equivariant. 

Lemma 4.2. The composition of a vertical element with an arbitrary element of Q x Q is G-equivariant, 

ihqaihgh^ 1 ) ■ {hq ,hq{) = h ■ i qo (g) -{q ,qi). 



Proof. Consider the following computation, 

ihqaihgh^ 1 ) ■ (hqv,hqi) 



= (hq ,hgh~ 1 hq 1 ) 
= (hqo,hgq!) 
= h(q ,gqi) 

= h ■ i qo (g) ■ (q ,qi) ■ □ 



Having made sense of how to compose an arbitrary element of Q x Q with a vertical element, we are in a 
position to introduce the notion of a discrete connection. 

A discrete connection is a G-equivariant choice of a subset of Q x Q called the discrete horizontal 
space, that is complementary to the discrete vertical space. In particular, given (qo,qi) G Q x Q, a discrete 
connection decomposes this into the horizontal component, hor(q , q\), and the vertical component, 
ver(go, qi), such that 

ver(g , qi) ■ hor^o, <7i) = (qo, Qi) , 
in the sense of the composition of a vertical element with an arbitrary element we defined previously. 
Furthermore, the G-equivariance condition states that 

hor(OTo,59i) = 9 ' hor(g , <Zi) , 

and 

verOft, gqi ) = 9 ■ ver(g , qi) ■ 

4.2. Discrete Atiyah Sequence. Recall that we obtain a short exact sequence corresponding to the de- 
composition of TQ into horizontal and vertical spaces. Due to the equivariant nature of the decomposition, 
quotienting this short exact sequence yields the Atiyah sequence. In this subsection, we will introduce the 
analogous discrete objects. 

Short Exact Sequence. The decomposition of the pair groupoid Q x Q, into discrete horizontal and 
vertical spaces, yields the following short exact sequence of bundles over Q. 

> VerQ^-^Q x Q -^-4 (7r,7r)*5 x S > 0, 

where VerQ is the discrete vertical subspace of Q x Q, and (ir, ir)*S x S is the pull-back of S x S by the 
projection (ir, ir) : Q x Q — > S x S. 

Discrete Atiyah Sequence. When the short exact sequence above is quotiented modulo G, we obtain an 
exact sequence of bundles over S, 

► G (Q x Q)/G SxS ► , 

which we call the discrete Atiyah sequence. Here, G is an associated bundle (see Definition 2.4). In 
particular, 

G = Q x G G= (Q x G)/G, 
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where the action of G on Q x G is given by g(q, h) — (gq, ghg^ 1 ), which is the natural discrete analogue of 
the adjoint action of g on Q x g. Furthermore, ttq : G — > S is given by 7Tg([<7, s]g) = ^(q)- 

The maps in the discrete Atiyah sequence i : G — > (Q x Q)/G, and (77, 7r) : (Q x Q)/G] — > 5 x 5, are given 

by 

#]g) = [?,5?]g = [«?(s)]g, 

and 

(7r,7r)([q ,gi] 9 ) = (7rqo,7rgi) • 



4.3. Equivalent Representations of a Discrete Connection. In addition to the discrete connection 
which arises from a G-equivariant decomposition of the pair groupoid Q x Q into a discrete horizontal and 
vertical space, we have equivalent representations in terms of splittings of the discrete Atiyah sequence, as 
well as maps on the unreduced short exact sequence. 

Maps on the Unreduced Short Exact Sequence. These correspond to discrete analogues of the 
connection 1-form, and the horizontal lift. 

• Discrete connection 1-form, Ad ■ Q x Q — > G. 

• Discrete horizontal lift, (■,')J:SxQ-»QxQ. 



Maps That Yield a Splitting of the Discrete Atiyah Sequence. 

• (7Ti , ^4d) : (Q x Q)/G — > G, which is related to the discrete connection 1-form. 

• (•, -) h : S x S — > (Q x Q)/G, which is related to the discrete horizontal lift. 



Relating the Two Sets of Representations. These two sets of representations are related in the following 
way: 

• The maps on the unreduced short exact sequence are equivariant, and hence drop to the discrete 
Atiyah sequence, where they induce splittings of the short exact sequence. 

• The maps that yield splittings of the discrete Atiyah sequence can be extended equivariantly to 
recover the maps on the unreduced short exact sequence. 

Furthermore, standard results from homological algebra yield an equivalence between the two splittings of 
the discrete Atiyah sequence. 

In the rest of this section, we will also discuss in detail the method of moving between the various 
representations of the discrete connection. The organization of the rest of the section, and the subsections 
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in which we relate the various representations arc given in the following diagram. 



§4.4 

discrete connection 1-form 
A d : Q x Q -> G 



ij4.fi 



§4.6 

splitting (connection 1-form) 
(ir,A d ):(QxQ)/G^G 



§4.1 

discrete connection 
hor : Q x Q — ► Hor 9 
ver : Q x Q — > Ver g 




§4.5 




§4.5 

discrete horizontal lift 
(■,■)>> -.SxS^QxQ 



§4.7 



§4.7 

splitting (horizontal lift) 
(■,-) h :SxS^(QxQ)/G 



4.4. Discrete Connection 1-Form. Given a discrete connection on a principal fiber bundle ir : Q — > Q/G, 
we can represent this as a Lie group- valued discrete connection 1-form, Ad ■ Q x Q — > G, which is a 
natural generalization of the Lie algebra-valued connection 1-form on tangent bundles, A : TQ — > g, to the 
discrete context. 

Discrete Connection 1-Forms from Discrete Connections. The discrete connection 1-form is con- 
structed as follows. Given an element of the Lie group g E G, the discrete generator map g i q (g) yields an 
isomorphism between G and Ver 9 for each q € Q. For each (q , qi) E Q x Q, we define Ad(qo, <Zi) to be the 
unique g E G such that i q (g) is equal to the vertical component of (</o>?i)- ln particular, this is equivalent 
to the condition that the following statement holds, 

(<7o,<7i) =i qo (A d (qo,qi)) ■ hor(g , t?i). 

Remark 4.1. It follows from the above identity that the discrete horizontal space can also be expressed as 

Hor go = {(Qo,Qi) EQ xQ \ hor(g ,gi) = (<7o,<7i)} 

= {(QO, q\) E Q X Q | A d (qo,qi) = e} . 

We will now establish a few properties of the discrete connection 1-form. 

Proposition 4.3. The discrete connection 1-form, Ad ■ Q x Q — > G, satisfies the following properties. 

(1) The 1-form is G-equivariant, that is, 

Ad° Lg = Ig O Ad, 

which is the discrete analogue of the G-equivariance of the continuous connection, AoTL g = Ad g oA. 

(2) The 1-form induces a splitting of the Discrete Atiyah sequence, that is, 

Ad(i q (g)) = Ad(qa,gq ) = g, 
which is the discrete analogue of A(£q) — £. 

Proof. The proof relies on the properties of a discrete connection, and the definition of the discrete connection 
1-form. 
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(1) The discrete connection 1-form satisfies the condition 

(q ,qi) = i qo {Ad(qo,qi)) ■ hor(g ,gi) • 
If we denote hor(qo, qi) by (go, q~i), we have that 

(q ,qi) = (q a ,A d (qo,qi)qi) ■ 

Similarly, we have, 

(gqo,gqi) = ig qo (A d {gqo,gqi)) ■ hor(g<? ,.ggi) 

= igq (A d (gq ,gqi)) ■ g ■ hor(g ,?i) 
= (e, -A d (gqo, gqi))(gqo, gqi) 

= (gqo,A d (gqo,gqi)gqi) , 

where we have used the G-equivariance of the discrete horizontal space. By looking at the expressions 
for gqi and q\, we conclude that 

A d (gq , gqijgqi = gqi 

= gA d {qo,qi)qi , 
A d (gqn, gqi)g = gA d {qo,qi) , 
A d (gqo,gqi) = gA d {qo,qi)g~ x , 

which is precisely the statement that Ad° L g = I g o Ad, that is to say that Ad is G-equivariant. 

(2) Recall that i q {g) is an element of the discrete vertical space. Since the discrete horizontal space 
is complementary to the discrete vertical space, it follows that ver(i q (g)) = i q {g). Then, by the 
construction of the discrete connection 1-form, Ad(i q (g)) is the unique element of G such that 

i q (A d (i q (g))) = ver(i q (g)) = i q (g) . 

Since i q is an isomorphism between G and the discrete vertical space, we conclude that Ad(i q (g)) = .9, 
as desired. □ 

The second result is equivalent to the map recovering the discrete Euler-Poincare connection when re- 
stricted to a G-fiber, that is, Ad(x, go, x, gi) — g\g^ ■ In particular, it follows that the map is trivial when 
restricted to the diagonal space, that is, A d (q,q) = e. 

The properties of a discrete connection are discrete analogues of the properties of a continuous connection 
in the sense that if a discrete connection has a given property, the corresponding continuous connection 
which is induced in the infinitesimal limit has the analogous continuous property. The precise sense in which 
a discrete connection induces a continuous connection will be discussed in §5.2. 

Discrete Connections from Discrete Connection 1-Forms. Having shown how to obtain a discrete 
connection 1-form from a discrete connection, let us consider the converse case of obtaining a discrete 
connection from a discrete connection 1-form with the properties above. We do this by constructing the 
discrete horizontal and vertical components as follows. 

Definition 4.4. Given a discrete connection 1-form, Ad ■ Q x Q — > G that is G-equivariant and induces a 
splitting of the discrete Atiyah sequence, we define the horizontal component to be 

hor(<? ,<7i) = iqodAdiqo^i))^ 1 ) ■ (q ,qi) ■ 
The vertical component is given by 

ver(<7o,<7i) = i qo {A d {qo, <7i)) ■ 

Proposition 4.4. The discrete connection we obtain from a discrete connection 1-form has the following 
properties. 
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(1) The discrete connection yields a horizontal and vertical decomposition of Q x Q, in the sense that 

(<10,qi) = vcr(g , qi) ■ hor(q , <?i) , 

for all (q ,qi) G Q x Q. 

(2) The discrete connection is G-equivariant, in the sense that 

hor(firg , 9Qi) = 9 ■ hor(q ,qi) , 



vev{gq , gqi) = g ■ vei(q , q{) . 

Proof. The proof relies on the properties of the discrete connection 1-form, and the definitions of the discrete 
horizontal and vertical spaces. 

(1) Consider the following computation, 

ver(<? ,<7i) ■ hor(q> ,q<i) = i qo {A d (q , qi)) ■ i qo ((A d (qo, qi))^ 1 ) ■ (?0,9i) 
= iqoiAdiqo^^iAdiqo^qi))^ 1 ) ■ (?0,9i) 

= iq ( e ) ■ (qa,qi) 
= (qo,qi) , 

where we used that i q is a homomorphism (see Lemma 4.1). 

(2) Wc compute, 

hor(.ggo,.Wi) = i gqa ((A d (gq , .99i)) _1 ) ' (59o,OTi) 

= hqa (g(A d (q , qi)r 1 g~ 1 ) ■ (gq ,gqi) 
= {e,g{Ad(qo,g.i))~ 1 g~ 1 )(gqo,9qi) 

= (gqo,g(A d (q ,qiyr 1 qi) 

= 9 ■ i q ((A d (qo,qi)y 1 ) ■ (<?o,9l) 
= g ■ hor(q ,<7i) , 

where we have used the fact that the composition of a vertical element with an arbitrary element is 
G-equivariant (see Lemma 4.2). Similarly, we compute, 

ver(g<7o,59i) = igq (A d (gq , gqi)) 

= iggoigAciiq^q^g- 1 ) 

= {gqo, gA d (qo, q^g^ 1 gq ) 

= g ■ (qo,A d (qo,qi)q ) 

= g ■ i go (A d (qo,qi)) 

= .9-ver(q ,gi). □ 



Local Representation of the Discrete Connection 1-Form. Since the discrete connection 1-form can 
be thought of as comparing group fiber quantities at different base points, we obtain the natural identity 
that 

A d (gq ,hqi) = M d (<Zo, 9i).9 _1 • 
In a local trivialization, this corresponds to 

Ad(x ,go,xi,gi) = giA d (x ,e,x 1 ,e)go 1 . 
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We define 

A(x , x\) = Ad(x ,e, xi,e), 
which yields the local representation of the discrete connection 1-form. 

Definition 4.5. Given a discrete connection 1-form, Ad ■ Q x Q — > G, its local representation is given 
by 

A d (x ,g ,xi,gi) = g^x^^x^g^ 1 , 

where 

A(x ,xi) = A d {x , e, xi,e). 
Lemma 4.5. The local representation of a discrete connection is G-equivariant. 
Proof. Consider the following computation, 

Ad(g(x a ,go),g(x 1 ,g 1 )) = A d ((x ,ggo), (x 1 ,gg 1 )) 
= ggiA(x Q , x 1 ){gg y 1 
= 5(giA(x ,a;i).g ( 7 1 ) 3 " 1 

= gA d ((x ,go), (a:i,5i))s _1 , 
which shows that the local representation is G-equivariant, as expected. □ 

Notice also that in the pure group case, where Q = G, this recovers the discrete Euler-Poincare connection, 
as we would expect, since the shape space is trivial. In particular, xq = x\ = e, which implies that 
A(x ,xi) = Ad{e,e,e,e) = e, and A d (g ,gi) = gig^ 1 ■ 

Example 4.1. As an example, we construct the natural discrete analogue of the mechanical connection, 
A : TQ — ► q, by the following procedure, which yields a discrete connection 1-form, Ad ■ Q x Q — > G. 

(1) Given the point (qo,qi) E Q x Q, we construct the geodesic path qoi : [0, 1] — > Q with respect to the 
kinetic energy metric, such that goi(0) = qo, and qoi(l) = qi. 

(2) Project the geodesic path to the shape space, xoi(t) = irqoi(t), to obtain the curve xoi on S. 

(3) Taking the horizontal lift of x m to Q using the connection A yields q m . 

(4) There is a unique g £ G such that goi(l) = g • <Zoi(l)- 

(5) Define A d (qo, Qi) = 9- 

This discrete connection is consistent with the classical notion of a connection in the limit that qi approaches 
qo, in the usual sense in which discrete mechanics on Q x Q converges to continuous Lagrangian mechanics 
on TQ. As mentioned before, this statement is made more precise in $5.2. 

4.5. Discrete Horizontal Lift. The discrete horizontal lift of an element (xo, x\) € S x S is the subset 
of Q x Q that are horizontal elements, and project to (xo,xi). Once we specify the base point q € Q, the 
discrete horizontal lift is unique, and we introduce the map (•,•)£ : S x S — > Q x Q. 

Discrete Horizontal Lifts from Discrete Connections. The discrete horizontal lift can be constructed 
once the discrete horizontal space is defined by a choice of discrete connection. 

Definition 4.6. The discrete horizontal lift is the unique map (■,■)£ : S x S —> Q x Q, such that 

(7r,7r) • (x ,xi)g = {x ,xi) , 

and 

{xo,Xi)g £ Hor 9 . 

Lemma 4.6. The discrete horizontal lift is G-equivariant, which is to say that 

(xo,xi)g q = g ■ (z ,2;i)£ . 
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Proof. Given (xq, xi) € S x S, denote (xo, x\ )^ by (qo, qi). Then, by the definition of the discrete horizontal 
lift, we have that 

(7r,7r) • (q ,qi) = {X ,X!) , 

and it follows that 

(7r,7r) • (gq ,gqi) = {x ,x 1 ). 
Also, from the definition of the discrete horizontal lift, 

(Qo,Qi) e Hor 90 , 
and by the G-equivariance of the horizontal space, 

(.990,. .991 ) G 9 ■ Hor^ = Hor g90 . 
This implies that (gqo,gqi) satisfies the conditions for being the discrete horizontal lift of (x ,xi) with 
basepoint gq . Therefore, (x , xi) h gqa = (.99o, 39i) = 9 • {Qo, Qi) = .9 ' {xo,xi)g , as desired. □ 

Discrete Connections from Discrete Horizontal Lifts. Conversely, given a discrete horizontal lift, we 
can recover a discrete connection. 

Definition 4.7. Given a discrete horizontal lift, we define the horizontal component to be 

hor(<7o,9i) = O(9o,9i))g > 
and the vertical component is given by 

ver(<?o,9i) = , 
where g is the unique group element such that 

(9o,9i) = ig (g) ■ hor(g ,9i) • 

The last expression simply states that the discrete horizontal and vertical space are complementary with 
respect to the composition we defined between a vertical element and an arbitrary element of Q x Q. 

Discrete Horizontal Lifts from Discrete Connection 1-Forms. We wish to construct a discrete 
horizontal lift (•, -) h : S x S — > (Q x Q)/G, given a discrete connection Ad : Q x Q — > G. We state the 
construction of such a discrete horizontal lift as a proposition. 

Proposition 4.7. Given a discrete connection 1-form, Ad ■ Q x Q — > G, the discrete horizontal lift is given 
by 

(x , Xl ) h = [7r _1 (a;o,a:i) n^ 1 (e)] G . 
Furthermore, the discrete horizontal lift satisfies the following identity, 

i qo {-Ad(qo,Qi)) ■ (7r(qo,9i))g = (9o,9i), 

which implies that the discrete connection 1-form and the discrete horizontal lift induces a horizontal and 
vertical decomposition of Q x Q. 

The horizontal lift can be expressed in a local trivialization, where q = (x ,go), using the local expression 
for the discrete connection, 

( x o,Xi)g = (xo,go,a;i,.9o(-4(a;o,a;i)) _1 )- 

Proof. We will show that this operation is well-defined on the quotient space. Using the local representation 
of the discrete connection in the local trivialization (see Definition 4.5), we have, 

A^ 1 (e)nTT^ 1 (x ,x 1 ) 

= {(xo,g,xi,g- (A(x ,xi)) _1 ) | x ,xi e S,g e G} 
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n {(xo,ho,xi,hi) | ho, hi e G} 

= {{xo,g,xi,g- (A(x ,xi))~ 1 ) \ h e G} 

= G ■ (x ,e,xi, (A(xo,a;i)) _1 ), 

which is a well-defined element of (Q x Q)/G. Since this is true in a local trivialization, and both the 
discrete connection and projection operators are globally defined, this inverse coset is globally well-defined 
as an element of (Q x Q)/G. 

In particular, the computation above allows us to obtain a local expression for the discrete horizontal lift 
in terms of the local representation of the discrete connection. That is, 

{xo,xi) h {xae) = (x a ,e,x 1 ,(A(x ,xi)y 1 ), 

(x ,xi) h = [(x ,e,a;i, (A(x , £i)) _1 )]g- 

By the properties of the discrete horizontal lift, this extends to 7r _1 (a;o, £1) cQxQ, 

{x ,xi) h (xQg) = {x ,xi) h g(xoe) 

= g ■ {xo,xi) h {xote) 

= g(x ,e,xi, (A(a;o,a;i)) _1 ) 

= {xo,g,xi,g(A(x ,Xi))~ 1 ). 

To prove the second claim, we have in the local trivialization of Q x Q, (qo,qi) = {xo,go,Xi,gi). Then, 
by the result above, 

W?o,gi))g = {n{qo,qi))( Xo ,g ) = {x , go, x u g (A(x a , xi))- 1 ). 
Also, by the local representation of the discrete connection, 

A d (q ,qi) = giA(xo,xi)gQ l . 

Therefore, 

i qo {Ad{qo,qi) ■ (n(qo,qi))q = (e,A d (qo,qi)) ■ {^(Qo,Qi))q 

= (e, giA(x , x^g^ 1 ) ■ (x , go,xi, g (A(x , x^y 1 ) 
= {xo,go,xi, (giA(x ,xi)g ~ 1 )(go(A(xo,xi))~ 1 )) 
= {xo,go,xi,gi) 
= (<7o,9i), 

as claimed. □ 



Discrete Connection 1-Forms from Discrete Horizontal Lifts. Given a horizontal lift (•, ')J :SxS^ 
Q x Q, we wish to construct a discrete connection 1-form, Ad ■ Q x Q — > G. 

Lemma 4.8. Given a discrete horizontal lift, (•,•)£ : S x S — »• Q x Q, the discrete connection 1-form, 

Ad ■ Q x Q — > G, is uniquely defined by the following identity, 

i qo {-A d (qo,qi)) ■ {n(qo,qi))q = (qo,qi)- 

Proof. To show that this construction is well-defined, we note that TTi(qo,qi) = 7r i( 7r (<Zo, 9i))q , by the 
construction of (•, -)g lrom ('>') h '■ S x S ^ (Q x Q)/G. Furthermore, ~K2{qo,qi) an d 7r 2 (7r((7o, qi))q are i n 
the same fiber of the principal bundle tt : Q — > Q/G and are therefore related by a unique element jeG. 
Since this element is unique, Ad{qo, qi) is uniquely defined by the identity. □ 
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4.6. Splitting of the Discrete Atiyah Sequence (Connection 1-Form). Consider the discrete Atiyah 
sequence, 

o ► g (Q x Q)/g Sxs ► o . 

Given a short exact sequence 

/ a 
> A 1 ^zr^ B M > , 

k h 

there are three equivalent conditions under which the exact sequence is split. They are as follows, 

(1) There is a homomorphism h : A<i — ► B with g oh — 1a 2 ] 

(2) There is a homomorphism k : B — > A\ with k o / = 

(3) The given sequence is isomorphic (with identity maps on A\ and A 2 ) to the direct sum short exact 
sequence, 

> A x Ai © A 2 — A 2 >■ , 

and in particular, B = A\ ffi A 2 . 
We will address all three conditions in this and the next two subsections. 

Splittings from Discrete Connection 1-Forms. A discrete connection 1-form, Ad ■ Q x Q — > G, induces 
a splitting of the discrete Atiyah sequence, in the sense that 

(iT!, A d ) oi = l 6 . 

Lemma 4.9. Given a discrete connection 1-form, Ad ■ Q x Q — > G, we obtain a splitting of the discrete 
Atiyah sequence, ip : (Q x Q)/G — > G, which is given by 

</?([<7o,<7i]g) = [qo,A d (qo,qi)]G- 

We denote this map by (ni,Ad)- 

Proof. This expression is well-defined, as the following computation shows, 

<p{\gqo, gqite) = [gqo,A d (gqo,gqi)]G 
= [gqo, gA d (qo, q^g^h 
= <p([qo,qi]G)- 

Furthermore, since 

(TT 1 ,A d ) oi([q,g} G ) = {TT 1 ,Ad){[q,gq] G ) 

= [ni(q,gq),Ad{q,gq)]G 

= [q,g]c, 

it follows that we obtain a splitting of the discrete Atiyah sequence. □ 



Discrete Connection 1-Forms from Splittings. Given a splitting of the discrete Atiyah sequence, we 
can obtain a discrete connection 1-form using the following construction. 

Given [<7o,9i]g G (Q x Q)/G, we obtain from the splitting of the discrete Atiyah sequence an element, 
[q-,g]G £ G. Viewing [q,g]c as a subset of Q x G, consider the unique g such that (qo,g) S [q,g]c C Q x G. 
Then, we define 

A d (x ,e,x l ,gQ 1 g 1 ) = g. 
We extend this definition to the whole of Q x Q by equivariance, 

■Ad(x ,go,xi,gi) ^goggo 1 - 
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Lemma 4.10. Given a splitting of the discrete Atiyah sequence, the construction above yields a discrete 
connection 1-form with the requisite properties. 

Proof. To show that the Ad satisfies the properties of a discrete connection 1-form, we first note that 
equivariance follows from the construction. 

Since we have a splitting, it follows that ip([q, gq]a) — [q,9]g, a s f composed with the map from G to 
(Q x Q)/G is the identity on G. Using a local trivialization, we have, 

[qo j 9} g = <p([qo,gqo]G) 

= ^([(x ,e),(x ,5(7 1 55o)]G) 

= [{xo,e),g] G 

= [{xo,go),go99o 1 ]G- 

Then, by definition, 

A d ((x , e), (x , g^ggo)) = g, 
and furthermore, g = 5o55o~ 1 ■ From this, we conclude that 

A d (qo,gqo) = A d ((x ,g ), (x ,55o)) 

= g A d ((x ,e), (xo^^ggo^go 1 

= 5o55(7 1 
= 5- 

Therefore, we have that Ad(qo,gqo) = g, which together with equivariance implies that Ad is a discrete 
connection 1-form. □ 

4.7. Splitting of the Discrete Atiyah Sequence (Horizontal Lift). As was the case with the discrete 
connection 1-form, the discrete horizontal lift is in one-to-one correspondence with splittings of the discrete 
Atiyah sequence, and they are related by taking the quotient, or extending by G-equi variance, as appropriate. 
Splittings from Discrete Horizontal Lifts. Given a discrete horizontal lift, we obtain a splitting by 

taking its quotient. 

Lemma 4.11. Given a discrete horizontal lift, (•, •)£ : S x £ — > Q x Q, the map (•, ■) fl :SxS-»(Qx Q)/G, 
which is given by 

induces a splitting of the discrete Atiyah sequence. 
Proof. We compute, 

(tt.tt) o (x ,a;i)' 1 = (7r,7r)([(x ,a; 1 )( i a;oie) ] G ) 
= {xo,xi) , 

where we used the G-equivariance of the discrete horizontal lift, and the property that (7r,7r) • (x ,xi)^ = 
(x ,xi) for any q e Q. This implies that {n,ir) o (•, -) h = lsxs, as desired. □ 

Discrete Horizontal Lifts from Splittings. Given a splitting, (•, -) h : S x S — ► (Q x Q)/G, we obtain a 
discrete horizontal lift, (•, •)£ : S x S — > Q x Q, using the following construction. 

We denote by (xo,Xi)g the unique element in (xo,xi) h , thought of as a subset of Q x Q, such that the 
first component is qo. This is the discrete horizontal lift of the point (xq, x\) € S x S where the base point 
is specified. 
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Lemma 4.12. Given a splitting of the discrete Atiyah sequence, the construction above yields a discrete 
horizontal lift with the requisite properties. 

Proof. Since the quotient space (Q x Q)/G is obtained by the diagonal action of G on Q x Q, it follows 
that if (xo,xi)g G (xo,xi) h C Q x Q, then g ■ (x ,xi)^ o G (x n ,xi) h C Q x Q. Since the first component of 
G ■ (x ,xi)q is gq , and g ■ (x a ,xi)^ o G (x ,xi) h C Q x Q, we have that 

(zo,a;i)gg = S-(ar ,a;i)£ , 

which is to say that the discrete horizontal lift constructed above is G-equi variant. 

Since (•, -) h is a splitting of the discrete Atiyah sequence, we have that (it, ir) o (-,-) h = IsxS, and this 

implies that any element in (x ,xi) h , viewed as a subset of Q x Q, projects to (xq 7 xi). Therefore, the 

discrete horizontal lift we constructed above has the requisite properties. □ 

4.8. Isomorphism between (Q x Q)/G and (5* x S) ©G. The notion of a discrete connection is motivated 
by the desire to construct a global diffeomorphism between (Q x Q)/G — > S and (S x S) © G — > S. This 
is the discrete analogue of the identification between TQ/G — > Q/G and T{Q/G) © § — > Q/G which is the 
context for Lagrangian Reduction in Cendra et al. [2001]. Since a choice of discrete connection corresponds 
to a choice of splitting of the discrete Atiyah sequence, we have the following commutative diagram, where 
each row is a short exact sequence. 







(<?,3<?) ( 7r , 7r ) 

-+ g (Q x Q)/g Sxs- 



■^0 



(ni, Ad) 







-> G G®(Sx5) 5xS^ 



-i-0 



Here, wc see how the identification between (Q x Q) /G and (S* x S) © G are naturally related to the discrete 
connection and the discrete horizontal lift. 

Recall that the discrete adjoint bundle G is the associated bundle one obtains when M = G, and p g acts 
by conjugation. Furthermore, the action of G on Q x Q is by the diagonal action, and the action of G x G 
on Q x Q is component- wise. 

Proposition 4.13. The map a^ d : (Q x Q)/G — > (S x S) © G defined by 

dA d {[qo,q\]G) = (^qo^qi) © [<?o,-4d(9o,?i)]G, 

is a well-defined bundle isomorphism. The inverse of ctA d is given by 

a^((a;o,a;i) © [q,g]o) = [(e,g) ■ (x , x^g, 
for any q G Q such that nq = x ■ 

Proof. To show that is well-defined, note that for any g e G, we have that 

{-xgqo^gqi) = (7rg ,7rgi), 

and also, 

[gqo,Ad{gqo,gqi)h = [gqo, gA d (q , qijg^h 
= [qo,A d (qo,qi)]G- 

Then, we see that 

aA d ([gqo,gqih) = a Ad {[Qo,Qi]G)- 

To show that a A ^ is well-defined, note that for any k G G, 
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and that 

a^ d ((x , x{) © [kq, kgk^] G ) = [(e, kgk" 1 ) ■ (x , Zi^Jg 

= [(e, kgk^ 1 ) ■ k ■ {x ,Xi)%] G 
= [(ek^kgk^k) ■ {x ,xi)^] G 
= [{ke,kg) ■ (x ,xi)g] G 
= [k ■ (e,g) ■ {x ,x{)\] G 
= [(e,g) ■ {x q ,xi)\]g 

= <x£ d ((x ,x 1 )®[q,g\ G )- □ 

Example 4.2. It is illustrative to consider the notion of a discrete connection, and the isomorphism, in 
the degenerate case when Q — G, which is the context of discrete Euler-Poincare reduction. Here, the 
isomorphism is between (G x G)/G and G, and the connection Ad ■ G x G — > G is given by 

•Ad{go,gi) =gi ■ go 1 - 

Then, we have that 

ou d ([0o, 9i]g) = (^0,^31) © [go,A d (qa,qi)]G 
= (e,e) © [ffo,5i5o • 

Taking the inverse, we have, 

auiGscSiffo^c) = [{e^igo 1 } • (e,e)Jj G 
= [(e,5i£fa • (9o,9o)]g 
= [ego,5iffo" 1 5o]G 
= \so,9i]g , 

as expected. 

4.9. Discrete Horizontal and Vertical Subspaces Revisited. Having now fully introduces all the equiv- 
alent representations of a discrete connection, we can revisit the notion of discrete horizontal and vertical 
subspaces in light of the new structures we have introduced. 
Consider the following split exact sequence, 

/ v g 
> Ai ^rr± B M > • 

k h 

We can decompose any element in B into a A\ and A 2 term by considering the following isomorphism, 

B = fok(B) ®hog(B). 

Similarly, in the discrete Atiyah sequence, we can decompose an element of (Q x Q)/G into a horizontal and 
vertical piece by performing the analogous construction on the split exact sequence 

o ► g (Q x Q)/g Sxs ► o . 

(Tl.^d) (-,-) h 

This allows us to define horizontal and vertical spaces associated with the pair groupoid Q x Q, in terms of 
all the structures we have introduced. 

Definition 4.8. The horizontal space is given by 

Uor q = {(q,q r ) eQxQ \ A d (q,q') = e} 
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= {(7r<7,xi)£ eQxQ\ Xl eS}. 
This is the discrete analogue of the statement Hor 9 = {v q G TQ \ A(v q ) = 0} = {(^7r 9 )g G TQ \ v^ q G TS}. 
Definition 4.9. The vertical space is given by 

Ver q = {{q,q') G Q x Q {n,w)(q,q') = e SxS } 
= {i g (g) | g G G}. 

This is the discrete analogue of the statement Ver q — {v q G TQ \ ir*(v q ) = 0} = {£g | £ G fl}. 

In particular, we can decompose an element of Q x Q into a horizontal and vertical component. 
Definition 4.10. The horizontal component o/(<7 ,9i) G Q x Q is given by 

hor(<7o,<7i) = ((■r) h ° K tt))(?o, ?i) = 0<?o, ^l)^- 
Definition 4.11. T7ie vertical component of (<?o,<Zi) G Q x Q is given 6t/ 

ver(q ,gi) = (i ° (tti, ^ d ))(g , 9i) = (qo,A d (q ,qi)q ) = i go (A d (q , qi))- 
Lemma 4.14. The horizontal component can be expressed as 

hor(<7o,<7i) = igodAdiqo^i))^ 1 ) ■ {qo,qi)- 

Proof. 

iqaiAdiqa^iY 1 ) ■ (qo,qi) = (qo, (A d {q , gi)) _1 <?o) • (qo,qi) 

= (e, (A d (qo, qi)y 1 )(qo,qi) 

= (q a , (Adiqo^iyy^l)- 
Clearly, (tt, n)(q , (Ad(qo, = (tt, 7r)(g , Furthermore, 

A d (q , (A d (qo,qi)y 1 qi) = (Ad(qo,qi)y 1 A d (qo,qi) = e. 
Therefore, by definition, (q , (A d (qo,qi)y 1 qi) = hor(g , qi)- □ 

Lemma 4.15. The horizontal and vertical operators satisfy the following identity, 

ver(<7o,<7i) • hor(g , <7i) = (qo,qi)- 

Proof. 

ver(q ,qi) ■ hor(g ,gi) = i qo (A d {q , qi)) ■ (i qo ((A d {q , qi))^ 1 ) ■ {qo,qi)) 
= {e,A d {q ,qi)){e, (A d (qo, 9i)) _1 )(9o, qi) 
= {e,A d {q ,qi)){q , (A d (qo, ?i))~V) 

= (?o, AK9o,gi)(Ai(<7o,<7i)rV) 
= (qo,qi), 

as desired. □ 

5. Geometric Structures Derived from the Discrete Connection 

In this section, we will introduce some of the additional geometric structures that can be derived from a 
choice of discrete connection. These structures include an extension of the pair groupoid composition to take 
into account the principal bundle structure, continuous connections that are a limit of a discrete connection, 
and higher-order connection-like structures. 
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5.1. Extending the Pair Groupoid Composition. Recall that the composition of a vertical element 
(qo,gqo) with an element (50,51) is given by 

(50,550) • (50,51) = (50,551)- 

The choice of a discrete connection allows us to further extend the composition, in a manner that is relevant 
in describing the curvature of a discrete connection. The decomposition of an element of Q x Q into a 
horizontal and vertical piece naturally suggests a generalization of the composition operation on Q x Q 
(viewed as a pair groupoid), by using the discrete connection, and the principal bundle structure of Q. 

We wish to define a composition on Q x Q such that the composition of (qo,qi) ■ (90, 5i) is defined 
whenever nqi = nq . Furthermore, we require that the extended composition be consistent with the vertical 
composition we introduced previously, as well as the pair groupoid composition, whenever their domains of 
definition coincide. 

The extended composition is obtained by left translating (q , qi) by a group element h, such that q\ = hqo, 
and then using the pair groupoid composition on (50,91) and the left translated term h(q~o,9i)- This yields 
the following intrinsic definition of the extended composition. 

Definition 5.1. The extended pair groupoid composition of (q , q{), (go, 9i) G Q x Q is defined whenever 
■nqi = nq , and it is given by 

(90,91) ■ (50,51) = {qo,-Ad(qo,qi)qi)- 

As the following lemmas show, this extended composition is consistent with the vertical composition and 
the pair groupoid composition. 

Lemma 5.1. The extended pair groupoid composition is consistent with the composition of a vertical element 
with an arbitrary element. 

Proof. Consider the composition of a vertical element with an arbitrary element, 

(qo,gqo) ■ (90,51) = (50,591)- 

This is consistent with the result using the extended composition, 

(qo,gqo) ■ (qo,qi) = {qo,A d (q ,gq )qi) 
= (50,591) , 

where we used that the discrete connection yields a splitting of the Atiyah sequence. □ 
Lemma 5.2. The extended pair groupoid composition is consistent with the pair groupoid composition. 
Proof. The pair groupoid composition is given by 

(5o,5i) • (51,52) = (50,52) • 
This is consistent with the extended composition, 

(5o,5i) • (51,52) = (50,^(51,51)52) 
= (5o,e52) 

= (5o,52)- □ 

The extended composition is G-equi variant, and is well-defined on the quotient space, as the following 
lemma shows. 

Lemma 5.3. The composition ■ : (Q x Q) x (Q x Q) — > (Q x Q) is G-equivariant, that is, 

(.55o,35i) • (35o,g5i) = g ■ ((5o,5i) ' (5o, 5i))- 
Furthermore, the composition induces a well-defined quotient composition ■ : {{Q x Q) x {Q x Q))/G —> 
(Q x Q)/G. 
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Proof. Given g G G, we consider, 

(gqo,gqi) ■ (Mo, Mi) = {gqo,A d {gq ,gqi)gqi) 
= (gqo,gA d (q ,qi)g~ 1 gqi) 
= (gqo,gA d {qo,qi)qi) 
= g ■ (qo,A d (q ,qi)qi) 
= g ■ ((50, qi) ■ (qo,qi)) , 

where we used the equivariance of the discrete connection. It follows that the composition is equivariant. 
Furthermore, 

[(gqo,gqi) ■ {gqo,gqi)h = \{qo,qi) ■ (?o,?i)]g, 

which means that • : {{Q x Q) x (Q x Q))/G — » (Q x Q)/G is well-defined. □ 
Corollary 5.4. TTie composition of n-terms is G- equivariant. That is to say, 

(gqlgql) ■ {gqlgql) ■■■■■ W 1 .^™" 1 ) ■ {9Qo,gq?) 

= g ■ ((go 1 .?! 1 ) • (?o^i) • • • • • KT^r 1 ) • («?,«?))■ 

Proof. The result follows by induction on the previous lemma. □ 

We find that the extended composition we have constructed on the pair groupoid is associative. However, 
as we shall see in §7.3, composing pair groupoid elements about a loop in the shape space will not in general 
yield the identity element cq x q, and the defect represents the holonomy about the loop, which is related to 
curvature. This may yield the discrete analogue of the expression giving the geometric phase in terms of a 
loop integral (in shape space) of the curvature of the connection. 

Lemma 5.5. The composition ■ : (Q x Q) x (Q x Q) — > (Q x Q) is associative. That is, 

(((&«?) • (qlql)) ■ (q 2 ,qi) = (gg,g?) • ((gUi) ■ («g,9?))- 
Proof. Evaluating the left-hand side, we obtain 

H,q°i) ■ (qlql)) ■ {qUD = («8,Ai(«5,«?)«i) • (qlql) 

= {q%,A d {q%,A d (ql,q1)ql)q1) 
= {ql,A d {ql,q\)A d {ql,q\)ql), 

and the right-hand side is given by 

(q° ,q°i) ■ ((qlQi) ■ {qlql)) = {qlql) ■ ialMfiA)® 

= {qlA d {ql,ql)A d {qlq\)ql). 

Therefore, 

{{qlql) ■ {qlql)) ■ {qlql) = {qlql) ■ {{qUD ■ {qlql)) , 

and the extended groupoid composition is associative. □ 

5.2. Continuous Connections from Discrete Connections. Given a discrete G-valued connection 1- 
form, A d : Q x Q — > G, we associate to it a continuous g-valued connection 1-form, A : TQ — > g, by the 
following construction, 

A{[q(-)]) = [A d (q(0),q{-m 
where [•] denotes the equivalence class of curves associated with a tangent vector. 

This uses the intrinsic representation of the tangent bundle, which is obtained by identifying a tangent 
vector at a point on the manifold with the equivalence class of curves on the manifold going through the 
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point, such that the tangent to the curve at the point is given by the tangent vector, which was illustrated 
earlier in Figure 9 on page 10. 

More explicitly, given v q E TQ, we consider an associated curve q : [0, 1] — > Q, and construct the curve 
g : [0, 1] -> G, given by 

g(t)=A d (q(0),q(t)). 

Then, 

AM = 37 »(*)■ 
Ul t=0 

When computing the equations in discrete reduction theory, it is often necessary to consider horizontal 
and vertical variations, which we introduce below. 

Definition 5.2. We introduce the vertical variation of a point (qo,qi) 6 Q x Q. Given a curve q\ : 
[0, 1] — * Q, such that gf (0) = qi, the vertical variation is given by 

ver(g ,9i) 



x d 
vcr oq = — 

ale 



e=0 



Definition 5.3. We introduce the horizontal variation of a point (qo,qi) G Q x Q. Given a curve 
q\ : [0, 1] — > Q, suc/i i/iai (0) = qi, the horizontal variation is given by 



hordq = -7- 
de 



_hov(q ,qt)=f e 



(«(qo,ql)) h qo 

e=0 



5.3. Connection-Like Structures on Higher-Order Tangent Bundles. Given a continuous connec- 
tion, we can construct connection-like structures on higher-order tangent bundles. This construction is de- 
scribed in detail in Lemma 3.2.1 of Cendra et al. [2001]. In particular, given a connection 1-form, A : TQ — > g, 
we obtain a well-defined map, A k : T^Q — > kg. 

As we will see later, these connection-like structures on higher-order tangent bundles will provide an 
intrinsic method of characterizing the order of approximation of a continuous connection by a discrete 
connection. 

We will describe the discrete analogue of this construction. To begin, the discrete analogue of the fc-th 
order tangent bundle, T^Q, is k + 1 copies of Q, namely Q k+1 . Intermediate spaces between T^Q and 
Q k+1 arise in the general theory of multi-spaces, which is introduced in Olver [2001]. 

The discrete analogue of tangent lifts, and their higher-order analogues, are obtained by componentwise 
application of the map, since a tangent lift of a map is computed by applying the map to a representative 
curve, and taking its equivalence class. Therefore, given a map / : M — > N, we have the naturally induced 
map, 

T^f : M k+1 - N k+1 given by f(m , . . . ,m k ) = (f(m ), . . . , f(m k )). 

And in particular, the group action is lifted to the diagonal group action on the product space. 
The discrete connection can be extended to Q k+1 in the natural way, A\ : Q k+1 — > ®i=qG = kG, 

A k d (q ,...,q k ) = ®t=oA d (q h qi +1 ). 
Similarly, we can define the map from Q k+1 to the Whitney sum of k copies of the conjugate bundle G 

by 

Q k+1 ^kG by (qo,...,qk)^®iJ[qo,A d (qi,qi + i)]G- 
In a natural way, we have the following proposition. 

Proposition 5.6. The map 

a A k : Q k+1 — > (Q/G) k+1 X-q/g kG 

defined by 

a A k{q Q , ...,q k ) = (7T<?o, • • .,7rg fe ) x Q/G (B k ^[qo, A d (qi, qi+i)] G , 
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is a well-defined bundle isomorphism. The inverse of a A k is given by 



a A l ((x , ...,x k ) -Xq/g ©to fa>9i] G ) 

= [(e, 5o, 5190, • • • ,9k-i ■ --go)) ■ {%o, ■ 



, x k) qo ]G, 



where (xo, ■ ■ ■ , Xk)q — (<7o, ■ ■ ■ , Qk) is defined by the conditions: 

Qo = qo, 
nqi = xi, 
A d (qi,qi+i) = e. 

Remark 5.1. In a local trivialization, where qo = (h ,x n ), we have, 

qi+i = ((A(xi,xi +1 ) ■ ■ ■ A(x , xi))" 1 ^, xi +1 ) . 



While we saw in the previous section that a discrete connection induces a continuous connection in 
the limit, we are often concerned with constructing discrete connections that approximate a continuous 
connection to a given order of approximation. This section will address the question of what it means for a 
discrete connection to approximate a continuous connection to a given order, as well as introduce methods 
for constructing such discrete connections. 

6.1. Exact Discrete Connection. It is interesting from the point of view of computation to construct an 
exact discrete connection associated with a prescribed continuous connection, so that we can make sense 
of the statement that a given discrete connection is a fc-th order approximation of a continuous connection. 

Additional Structure. To construct the exact discrete connection, we require that the configuration mani- 
fold Q be endowed with a bi-invariant Riemannian metric, with the property that the associated exponential 
map, 

exp :TQ-,Q, 
is consistent with the group action, in the sense that 

exp(£g(gO) = cxp(£) • q. 
We extend the exponential to Q x Q as follows, 

exp :TQ^QxQ, 

v q >-> (q,exp(v g )), 

and denote the inverse by log : Q x Q — > TQ, which is defined in a neighborhood of the diagonal of Q x Q. 

Construction. Having introduced the appropriate structure on the configuration manifold, we define the 
exact discrete connection as follows. 

Definition 6.1. The Exact Discrete Connection associated with a prescribed continuous connection 
A : TQ — > g and a Riemannian metric is given by 



6. Computational Aspects 



A% (qo,qi) = cxp(^4(log(( ? o, <7i)))- 



This construction is more clearly illustrated in the following diagram, 



QxQ 




A 



exp 



!<>» 
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Properties. The exact discrete connection satisfies the properties of the discrete connection 1-form, in that 
it is G-equi variant, and it induces a splitting of the discrete Atiyah sequence. The equivariance of the exact 
discrete connection arises from the fact that each of the composed maps is equivariant, and the splitting 
condition, 

Af{i q {g))=g, 
is a consequence of the compatibility condition, 

exp(£<2(<?)) = cxp(£) • q. 

Since the logarithm map is only defined on a neighborhood of the diagonal of Q x Q, it follows that the 
exact discrete connection will have the same restriction on its domain of definition. 



Example 6.1 (Discrete Mechanical Connection). The continuous mechanical connection is defined by the 
following diagram, 



T*Q 

FL 

TQ- 



->9 



Correspondingly, the discrete mechanical connection is defined by the following diagram, 

£d 

r \ 

¥ Ld J ^ 

I 

>9 



Q x Q. 



cxp 



>G 

J 



A d 

This is consistent with our notion of an exact discrete mechanical connection as the following diagram 
illustrates, 

Jd 

\. 



Q x Q 



QxQ 



log 



T*Q 



FL I 



TQ 



->0 



cxp 



where the portion in the dotted box recovers the continuous mechanical connection. In checking G- equivariance, 
we use the equivariance o/exp : g — > G, Jd ■ Q x Q — > g* , and the equivariance ofl in the sense of a map 
I : Q -> L(g, g*), namely, l(gq) ■ Ad g £ = Ad*-i l(q) ■ £. 



6.2. Order of Approximation of a Connection. We have the necessary constructions to consider the 
order to which a discrete connection approximates a continuous connection. There are two equivalent ways 
of defining the order of approximation of a continuous connection by a discrete connection, the first is more 
analytical, and is given by the order of convergence in an appropriate norm on the group. 
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Definition 6.2 (Order of Connection, Analytic). A discrete connection Ad is a k-th order discrete connection 
if, k is the maximum integer for which 

30 < c < oo, 
3h Q > 0, 

such that 

sup \\A% (q,ex.p(hv q ))(A d (q,exp(hv q )))- 1 \\ < ch k+1 , Vh<h . 

v q e TQ, 
\v q \ = 1 

The second definition is more intrinsic, and is related to considering the infinitesimal limit of a discrete 
connection to connection- like structures on higher-order tangent bundles, without the need for the introduc- 
tion of the exact discrete connection. 

Recall from §5.2 that we can construct a continuous connection from a discrete connection by the following 
construction, 

A([q(-)]) = [A d (q(0), it))}- 
Given A k : Q k+1 — > kG, we can obtain the continuous limit A k : T^Q — > fcg in a similar fashion. 

Definition 6.3 (Order of Connection, Intrinsic). A discrete connection Ad is a k-th order approximation 
to A if, k is the maximum integer for which the diagram holds, 

A d :QxQ^G A : TQ — > g 

V v 

A\ : Q k+1 -^kG ^A k : T {k) Q -» kg 

Here, the double arrows represent the higher-order structures induced by the connections, and the dotted 
arrow represents convergence in the limit. 

6.3. Discrete Connections from Exponentiated Continuous Connections. To apply the exponential 
and logarithm approach to construct a discrete connection from a prescribed connection, in the sense of the 
diagram, 

Ad 

f \ 

we can rely on explicit expressions for the exponential and logarithm, or we can rely on approximations to 
the exponential and logarithm. 

The explicit formulas for the special Euclidean group are particularly useful for applying the theory of 
discrete connections to the geometric control of problems such as robotic manipulators, and clusters of 
satellites. In dealing with other configuration manifolds, approximants to the exponential and logarithm 
may be required due to the absence of explicit formulas. 

Even when explicit formulas are available, it may be desirable to rely on a more computationally efficient 
approximation, such as the Cayley transformation, methods based on Pade approximants (see, for exam- 
ple, Cardoso and Silva Leite [2001], Higham [2001]), or Lie group techniques (see, for example, Celledoni 
and Iserles [2000, 2001], Zanna and Munthe-Kaas [2001/02]). Clearly, these will yield a discrete connection 
that has an order of approximation equal to the lower of the two orders of approximation of the numerical 
schemes used for the exponential and the logarithm. 

When used in the context of geometric control, high-order approximations to the continuous connection 
may not be necessary, since the optimal trajectory is often recomputed at each step, and in such instances, 
a low-order approximation suffices. 
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6.4. Discrete Mechanical Connections and Discrete Lagrangians. We will introduce a discrete me- 
chanical connection that is consistent with the structure of discrete variational mechanics, and will yield a 
discrete connection that has an order of approximation that is equation to that obtained from the discrete 
mechanics. 

Consider a G-invariant fc-th order discrete Lagrangian, L d : Q x Q — > M, which is to say that 

L d (gqo,gqi) = L d (q ,qi), 

and 

L d = Lf^ + 0(h k+1 ), 
where the exact discrete Lagrangian, £° xact : Q x Q — > R, is given by 

,-h 



LT aCt (qo,qi) = I L(q 01 (t),q 01 (t))dt. 
Jo 



Here, qoi : [0, ft] — > Q is the solution of the Euler-Lagrange equations with (701 (0) = qo, and 901 (ft) = 3i- 
The exact discrete Lagrangian is a generator of the symplectic flow, coming from the Jacobi solution of the 
Hamilton-Jacobi equation. 

This fc-th order discrete Lagrangian yields a fc-th order accurate numerical update scheme, through the 
discrete Euler-Lagrange equations, 

D 2 L d (q ,qi) + D 1 L d {q 1 , q 2 ) = 0, 

which implicitly defines a discrete flow $ : (go, ?i) >-> (<Zi, q 2 )- By pushing this numerical scheme forward to 
T*Q using the discrete fiber derivative ¥L d : Q x Q — > T*Q, which maps (qo, qi) 1— > (g 0) —D\L d {qQ, qi)), we 
can obtain a Symplectic Partitioned Runge-Kutta scheme of the same order. 
We also introduce the discrete momentum map, J d : Q x Q — > g* , given by 

(Jd(qo,qi),Q = -£>iL d (go,9i) • CqM- 

The discrete Lagrangian is G-invariant, which implies that for any £ s 0, we have, 
L d (qo,qi) = Ld(exp(£i) • <7 ,cxp(£i) ■ gi), 



dt 



L d (exp(£i) • g ,exp(£i) • ft) 

t=o 

= DxL^qn, qi) ■ Cq(<7o) + D 2 L d (q Q , ft) ■ £q(ft). 
If we restrict to the flow of the discrete Euler-Lagrange equations, we have that 

{DiL d {q a , qi ) + D 2 L d q(q 1 ,q 2 )) -£q(<Zi) =0, 
which upon substitution into the previous equation, yields 

L>ii d (<7o,<7i) -£q(<Zo) - D 1 L d (q 1 ,q 2 ) -?q(<?i) = 0, 

-£>iL d (gi,g 2 ) -£g(ft) = --t>i£d(<7o,<7i) -£q(<7o), 

(Jd(qi,q2),tQ(qi)) = (Jci{qo,qi),^q(qo)), 

which is the statement of the discrete Noether theorem, that the discrete momentum map is preserved by 
the discrete Euler-Lagrange flow. 

We note that the mechanical connection corresponds to a choice of horizontal space corresponding to the 
zero momentum surface. That is to say that the horizontal distribution corresponding to the continuous 
mechanical connection is 

Hor 9 = {v q e TQ I J(v q ) = 0}, 
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and the discrete horizontal subspace corresponding to the discrete mechanical connection is 

Hor g = {(<?> g') zQ xQ \ Jd(q, q') = 0} 

Remark 6.1. For the discrete horizontal subspace we defined above to have the correct dimensionality, 
the discrete Lagrangian needs to satisfy certain non- degeneracy conditions, which dictates the size of the 
neighborhood of the diagonal that the discrete connection is defined on. 

Since the continuous momentum map is preserved by the continuous Euler-Lagrange flow, and the discrete 
momentum map is preserved by the discrete Euler-Lagrange flow, it follows that the order of approximation 
of the continuous mechanical connection by the discrete mechanical connection is equal to the order of 
approximation of the continuous Euler-Lagrange flow by the discrete Euler-Lagrange flow. To construct a 
discrete mechanical connection of a prescribed order, we use the following procedure. 

(1) Consider a fc-th order G-invariant discrete Lagrangian, Ld : Q x Q — > R, 

L d = Lf ct + 0(h k+1 ). 

(2) Construct the corresponding discrete momentum map, Jd(qo, <Zi) - ► 0*, 

(Jd(qo,qi),Q = -DiL d {qo,qi) ■ Cq(<?o)- 

(3) Then, the fc-th order discrete mechanical connection is given implicitly by considering the condition 
for the discrete horizontal space, 

Ad(qo,qi) = e iff J d (q Q ,qi) = 0, 

and then extending the construction to the domain of definition by G-equivariance. 

(4) More explicitly, given (qo,qi) £ Q x Q, we consider a local trivialization, in which (qo,qi) — 
(xo, .go, xi,gi), and we find the unique g £ G such that 

Jd(xo, go,xi,g) = 0. 

Then, we have that 

■Ad(x ,g ,xi,g) = e, 

from which we conclude that 

■Ad(xo,9o,Xi,9i) = ■A d (xo,go,x 1 ,g 1 g~ 1 g) 
= gig^ 1 ■ A d (x ,go,xi,g) 
= gig^ 1 ■ 

7. Applications 

This section will sketch some of the applications of the mathematical machinery of discrete connections 
and discrete exterior calculus to problems in computational geometric mechanics, geometric control theory, 
and discrete Riemannian geometry. 

7.1. Discrete Lagrangian Reduction. Lagrangian reduction, which is the Lagrangian analogue of Poisson 
reduction on the Hamiltonian side, is associated with the reduction of Hamilton's variational principle for 
systems with symmetry. 

The variation of the action integral associated with a variation in the curve can be expressed in terms 
of the Euler-Lagrange operator, £C : T^Q -> T*Q. When the La grangian is G-invariant, the associated 
Euler-Lagrange operator is G-equivariant, and this induces a reduced Euler-Lagrange operator, [££]g ■ 
T^Q/G — > T*Q/G. The choice of a connection allows us to construct intrinsic coordinates on T^/G and 
T*Q/G, and the representation of the reduced Euler-Lagrange operator in these coordinates correspond to 
the Lagrange-Poincare operator, CP : T^(Q/G) x Q/G 2§ -> T*(Q/G) © §*. 
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The reduced equations obtained by reduction tend to have non-canonical symplectic structures. As such, 
naively applying standard symplectic algorithms to reduced equations can have undesirable consequences for 
the longtime behavior of the simulation, since it preserves the canonical symplectic form on the reduced space, 
as opposed to the reduced (non-canonical) symplectic form that is invariant under the reduced dynamics. 

This sends a cautionary message, that it is important to understand the reduction of discrete variational 
mechanics, since applying standard numerical algorithms to the reduced equations obtained from continuous 
reduction theory may yield undesirable results, inasmuch as long-term stability is concerned. 

Discrete connections on principal bundles provide the appropriate geometric structure to construct a 
discrete analogue of Lagrangian reduction. We first introduce the discrete Euler-Lagrange operator, which 
is constructed as follows. 

Discrete Euler Lagrange Operator. The discrete Euler-Lagrange operator, £ C d ■ Q 3 — ► T*Q satisfies 
the following property, 

de d (L d ) ■ <5q = ^2££ d (L d )(q k _ 1 ,q k ,q k+1 ) ■ 8q k . 
In coordinates, the discrete Euler-Lagrange operator has the form 

[D 2 L d (q k - 1 ,q k ) + D 1 L d (q k ,q k+1 )} dq k . 

Discrete Lagrange Poincare Operator. The map £C d {L d ) : Q 3 — > T*Q, being G-equi variant, induces 
a quotient map 

[£C d (L d )] G : Q 3 /G - T*Q/G, 
which depends only on the reduced discrete Lagrangian l d : (Q x Q)/G — ► R. We can therefore identify 
[£ C d (L d )] G with an operator £C d {l d ) which we call the reduced discrete Eulei — Lagrange operator. 

If in addition to the principal bundle structure, we have a discrete principal connection as described in 
the previous section, we can identify 

Q 3 /G with (Q/G) 3 x Q/G (G®G). 

The isomorphism between these two spaces is a consequence of Proposition 5.6, which is higher-order general- 
ization of Proposition 4.13. The discrete mechanical connection which was introduced in §6.4 is a particularly 
natural choice of discrete connection, since it does not require any ad hoc choices, as it is constructed directly 
from the discrete Lagrangian. 

Furthermore, each discrete G-valued connection 1-form, A d : Q x Q — > G, induces in the infinitesimal 
limit a continuous g-valued connection 1-form, A : TQ — ► g, as shown in §5.2. This continuous principal 
connection allows us to identify 

T*Q/G with T*(Q/G)©§*. 

The discrete Lagrange-Poincare operator, CV d (l d ) : (Q/G) 3 X-q/g {G © G) — > T*(Q/G) © §*, is obtained 
from the reduced discrete Euler-Lagrange operator by making the identifications obtained from the discrete 
connection structure. 

The splitting of the range space of£V d (l d ) as a direct product (as in §3.3 of Cendra et al. [2001]) naturally 
induces a decomposition of the discrete Lagrange-Poincare operator, 

CV d (l d ) = Kov{£V d (l d )) © VeT(£.r d {ld)), 

and this allows the discrete reduced equations to be decomposed in horizontal and vertical equations. 



7.2. Geometric Control Theory and Formations. There are well-established control algorithms for 
actuating a control system to achieve a desired reference configuration. In many problems of practical 
interest, the actuation of the mechanical system decomposes into shape and group variables in a natural 
fashion. 
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A canonical example of this is a satellite in motion about the Earth, where the orientation of the satellite is 
controlled by internal rotors through the use of holonomy and geometric phases, and the position is controlled 
by chemical propulsion. 

In this example, the configuration space is SE(3), the group is SO(3), and the shape space is M 3 . The group 
variable corresponds to the orientation, and the shape variable corresponds to the position. When given an 
initial and desired configuration, it is desirable, while computing the control inputs, to decompose the relative 
motion into a shape component and a group component, so that they can be individually actuated. 

Since the discrete connection is used here to provide an efficient choice of local coordinates for optimal 
control, the discrete connection is most naturally obtained by exponentiating the continuous connection, in 
the manner described in §6.3, in conjunction with the explicit formulas for the exponential and logarithm 
for SE(3). The natural choice of the continuous connection is one in which the horizontal space is given by 
the momentum surface corresponding to the current value of the momentum. 

To illustrate why it may not be desirable from a control-theoretic point of view to decompose the space 
using a trivial connection, consider a satellite that is in a tidally locked orbit about the Earth, with the 
initial and desired configuration as illustrated in Figure 10. 




Initial configuration Desired configuration 

Figure 10: Application of discrete connections to control. 

Here, if we choose a trivial connection, then the relative group element would be a rotation by 7r/4, but 
this choice is undesirable, since the motion is tidally locked, and moving the center of mass to the new 
location would result in a shift in the orientation by precisely the desired amount. In this case, the optimal 
control input should therefore only actuate in the shape variables, and the relative group element assigned 
to this pair of configurations should be the identity clement. 

As such, the extension of mechanically relevant connections to pairs of points in the configuration space 
with finite separation, through the use of a discrete connection, can be of immense value in geometric control 
theory. 

Similarly, in the case of formations, discrete connections allow for the orientation coordination problem to 
be handled in a more efficient manner, by taking into account the dynamic coupling of the shape and group 
motions automatically through the use of the discrete mechanical connection. 

7.3. Discrete Levi-Civita Connection. Vector bundle connections can be cast in the language of con- 
nections on principal bundles by considering the frame bundle consisting of oriented orthonormal frames 
over the manifold M, which is a principal SO(n) bundle, as originally proposed by Cartan [1983, 2001]. For 
related work on discrete connections on triangulated manifolds with applications to algebraic topology and 
the computation of Chern classes, please see Novikov [2003]. 

To construct our model of a discrete Riemannian manifold, we first trivialize the frame bundle to yield 
SO(n) x M. Then, Q = SO(n) x M, and G = SO(n). 
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Here, we introduce the notion of a semidiscretizcd principal bundle, where the shape space, S = Q/G, is 
discrctized as a simplicial complex K, and the structure group G remains continuous. In this context, the 
semidiscretization of the trivialization of the frame bundle is given by Q — SO(n) x K. 

A discrete connection is a map Ad ■ Q x Q — > G, and we can construct a candidate for the Levi-Civita 
connection on a simplicial complex K, using the discrete analogue of the frame bundle described above. 
However, we will first introduce the notion of a discrete Riemannian manifold. 

Definition 7.1. A discrete dual Riemannian manifold is a simplicial complex where each n-simplex 
a n is endowed with a constant Riemannian metric tensor g, such that the restriction of the metric tensor to 
a common face with an adjacent n-simplex is consistent. 

This is referred to as a discrete dual Riemannian manifold as we can equivalently think of associating 
a Riemannian metric tensor to each dual vertex, and as we shall see, by adopting Cartan's method of 
orthogonal frames (see, for example, Cartan [1983, 2001]), the connection is a SO(n)-valued discrete dual 
1-form, and the curvature is a SO(n)-valued discrete dual 2-form. 

For each n-simplex o~ n , consider an invcrtiblc transformation / of R" such that f*g = I. In the orthonormal 
space, we have a normal operator that maps a (n — l)-dimcnsional subspace to a generator of the orthogonal 
complement, denoted by _L. Then, we obtain a normal operator on the faces of a n by making the following 
diagram commute. 




The coordinate axes in the diagram represent the normalized eigenvectors of the metric, scaled by their 
respective eigenvalues. 

The local representation of the discrete connection is given by 

A d (K,Ro),(°i,Ri)) = RiA(^,a[ l )R \ 

and so the discrete connection is uniquely defined if we specify A(aJ},ai), where aj} and a" are adjacent 
n-simplices. Since they are adjacent, they share a (n — l)-simplex, denoted a™ -1 . In particular, this can 
then be thought of as a SO(n)-valued discrete dual 1-form, since to each dual 1-cell, ★cr™ -1 , we associate an 
element of SO(n). 

This element of SO(n) is computed as follows. 

(1) In each of the n-simplices, we have a normal direction associated with a™ -1 , denoted by 

_L (CT n_1 ,<7?) e K". 

(2) If these two normal directions are parallel, we set 

(A,*a n - 1 )=I, 

otherwise, we continue. 
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(3) Construct the (n — 2)-dimcnsional hyperplane P n 2 , given by the orthogonal complement to the 
span of the two normal directions. 

P"- 2 =1 (span(± (<T n -\aZ),±(a n -\a?))). 

(4) If *cr™~ 1 is oriented from CTq to tr™, we set 

(^W"- 1 ) = {Re SO(n) | R\ pn -2=I pn -,,R(± (tr"" 1 ,^)) =_L K^O}. 

The curvature of this discrete Levi-Civita connection is then a SO(n)-valued discrete dual 2-form. There 
is however the curious property that the boundary operator for a dual cell complex may not necessarily agree 
with the standard notion of boundary, since that may not be expressible in terms of a chain in the dual cell 
complex. This is primarily an issue on the boundary of the simplicial complex, and if we are in the interior, 
this is not a problem. 

Since curvature is a dual 2-form, it is associated with the dual of a codimension-two simplex, given by 
*<7™~ 2 . Consider the example illustrated in Figure 11. 




Figure 11: Discrete curvature as a discrete dual 2-form. 

The curvature B of the discrete Levi-Civita connection is given by 

{B,*a n - 2 ) = (cL4, *cr"- 2 ) = (A,d*a n - 2 ). 

As can be seen from the geometric region *d * er™~ 2 , the curvature associated with *er™~ 2 is given by the 
ordered product of the connection associated with the dual cells, ★ct™ -1 , where cr™ -1 >- a n ~ 2 . This also 
suggests that we can also think of the discrete connection as a primal (n — l)-form, and the curvature as a 
primal (n — 2)-form, where the curvature is obtained from the connection using the codifferential. 

When the group G is nonabelian, we see that the curvature is only defined up to conjugation, since we 
need to specify a dual vertex on the dual one-chain d-ko n ~ 2 from which to start composing group elements. 
To make this well-defined, we can adopt the approach used in defining the simplicial cup product, and assume 
that there is a partial ordering on dual vertices, which would make the curvature unambiguous. 

As a quality measure for simplicial triangulations of a Ricmannian manifold, having the curvature defined 
up to conjugation may be sufficient if we have a norm on SO(n) which is invariant under conjugation. As 
an example, taking the logarithm to get an element of the Lie algebra so(n), and then using a norm on this 
vector space yields a conjugation-invariant norm on SO(n). This allows us to detect regions of the mesh 
with high curvature, and selectively subdivide the triangulation in such regions. 

Similarly, we can define a discrete primal Riemannian manifold, where the Riemannian metric tensor is 
associated with primal vertices, and the connection is a G-valued primal 1-form, and the curvature is a 
G-valued primal 2-form. 

Abstract Simplicial Complex with a Local Metric. In Desbrun et al. [2003], the notion of an abstract 
simplicial complex with a local metric, defined on pairs of vertices that are adjacent, was introduced. 
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In this situation, we can compute the curvature around a loop in the mesh using local embcddings. We 
start with an initial n-simplex, which we endow with an orthonormal frame. By locally embedding adjacent n- 
simplices into Euclidean space, and parallel transporting the orthonormal frame, we will eventually transport 
the frame back to the initial simplex. 

The relative orientation between the original frame and the transported frame yields the integral of the 
curvature of the surface which is bounded by the traversed curve. This results from a simple application of 
the Generalized Stokes' theorem, and the fact that the curvature is given by the exterior derivative of the 
connection 1-form. 

8. Conclusions and Future Work 

We have introduced a complete characterization of discrete connections, in terms of horizontal and vertical 
spaces, discrete connection 1-forms, horizontal lifts, and splittings of the discrete Atiyah sequence. Geometric 
structures that can be derived from a given discrete connection have been discussed, including continuous 
connections, an extended pair groupoid composition, and higher-order analogues of the discrete connection. 
In addition, we have explored computational issues, such as order of accuracy, and the construction of discrete 
connections from continuous connections. 

Applications to discrete reduction theory, geometric control theory, and discrete geometry, have also been 
discussed, and it would be desirable to systematically apply the machinery of discrete connections to these 
problems. 

In addition, connections play a crucial role in representing the nonholonomic constraint distribution in 
nonholonomic mechanics, particularly when considering nonholonomic mechanical systems with symmetry, 
wherein the nonholonomic connection enters (see, for example, Bloch [2003] ) . There has been recent progress 
on constructing nonholonomic integrators in the work of Cortes [2002] and McLachlan and Perlmutter [2003] , 
but an intrinsically discrete notion of a connection remains absent from their work, and they do not consider 
the role of symmetry reduction in discrete nonholonomic mechanics. 

It would be very interesting to apply the general theory of discrete connections on principal bundles 
to nonholonomic mechanical systems with symmetry, and to cast the notion of a discrete nonholonomic 
constraint distribution and the nonholonomic connection in the language of discrete connections, and thereby 
develop a discrete theory of nonholonomic mechanics with symmetry. This would be particularly important 
for the numerical implementation of geometric control algorithms. 

The role of discrete connections in the study of discrete geometric phases would also be an area worth 
pursuing. A discrete analogue of the rigid-body phase formula, that involves the discrete mechanical connec- 
tion, that is exact for rigid-body simulations that use discrete variational mechanics, would yield significant 
insights into the geometric structure-preservation properties of variational integrators. In particular, it would 
provide much needed insight into how discretization interacts with geometric phases, and yield an under- 
standing how much of the phase drift observed in a numerical simulation is due to the underlying geometry 
of the mechanical system, and how much is due to the process of discretizing the system. 
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